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^< Abstract 

(N 



This paper examines an averaging technique apphcd to the trans- 
port equations as an alternative to vanishing viscosity. Such techniques 
have been shown to be vahd shock-regularizations of the Burgers equa- 
tion and the Euler equations, but has yet to be apphed to the similar 
transport equations. However, for this system, the classical notion of 
weak solutions is not always sufficient thus a more general notion of 
a distribution solution containing Dirac-delta functions must be intro- 



r^ duced. Moreover, the distribution solution to the Riemann problem is 

known to be the weak-* limit of the viscous perturbed transport equa- 
tions as viscosity vanishes. In comparison to the classical method of 
vanishing viscosity, the Riemann problem is examined for the averaged 
transport equation and it is shown that the same delta-shock distribu- 
tion solution is captured as filtering vanishes. Both mathematical and 
physical motivation are provided for the averaging techniques consid- 
ered including an existence and uniqueness result of solutions within 
the class of smooth initial conditions. 



1^ 

Q 1 Introduction 

O 

^ Consider the transport (or transportation) equations in multiple dimensions 

'k^ in which the continuity equation is adjoined with the inviscid Burgers equa- 

;_] tion 

^ f Pt + V • (pu) = 

Ui-Fu-Vu = (1.1) 

(p(x,0),u(a;,0)) = {po{x),uo{x)). 

The transport equations provide a simple system of two conservation laws 
for which no classical weak solutions may exist for some initial data, so a 
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distribution solution consisting of Dirac delta functions or a delta-shock so- 
lution must be introduced [HE]- On the other hand, the transport equations 
model the dynamics of particles that adhere to one another upon collision 
and has been studied as a simple cosmological model for describing the non- 
linear formation of large-scale structures in the Universe [3|. For example, 
u may represent the flow field carrying dust particles with density p, and 
the delta-shock wave represents a concentration of dust on a shock which 
attract the dust [5]. For the Riemann problem, the classical regularization 
of adding dissipation such as the adhesion model [3l [5] , 

Pt + V ■ (pu) = .^2) 

ut + u • Vu = i/Au, ^ ■ ^ 

have been shown to recover the delta-shock solution as i/ — )• 0. Similarly, 
K.T. Joseph |6j considered the Riemann problem for the following regular- 
ization 

Pt + {pu)x = eu^x ., „. 

ut + (-2-)x = e-Uxx 

and proved that the limit of the corresponding solutions (p*^, u^) may contain 
delta- measures as e — )■ 0. The authors in [2] examined a similar system in 
one dimension 

pt + {pu)x = , . 

ut + {u% = ^'-^^ 

in order to mathematically justify delta-shock waves. They, too, introduce 
viscous perturbations. One case involves adding the typical Laplacian term 

Pt + {pn)x = 
Ut + (n^jx = euxx 

while the other introduces time-dependent viscosity 

Pt + {pu)x = 

Ut + [u^)x = etuxx- 

They construct solutions to these perturbed systems and show that the 
delta-shock waves are the weak-* limits of the the solutions to the viscous 
systems as viscosity vanishes. 

The principle focus of this paper is to introduce an alternative regulariza- 
tion of the transport equations in which the delta-shock solution is captured 
similar to what vanishing viscosity achieves, yet without largely altering the 
initial structure of the equations. This particular regularization will utilize a 
spatial averaging of the nonlinear terms by convoluting them by some aver- 
aging kernel or low-pass filter rather than adding viscous terms. The notion 



of averaging the transport equations directly stems from the successful reg- 
ularization of the Burgers equation via averaging of the convective velocity 

m, 

Ut + UUx = 

u = g'^ * u (1-7) 

Their numerical study suggests that this filtering of the nonlinear term at- 
tenuates shock formation without adding viscosity or dispersion. Mohseni 
and Norgard [SI [9] as well as Bhat and Fetecau |10j later studied this regu- 
larization in one and multiple dimensions. It is shown that spatial averaging 
of the convective velocity properly regularizes the Burgers equation while 
numerical simulations suggest the entropy solution is captured as the filter- 
ing parameter, a, vanishes. In fact, the entropy solution is captured for a 
class of 'bell-shaped' C^ initial conditions [9]. 

It should be noted that this was not the first implementation of filtered con- 
vective velocities. In 1934, Jean Leray proposed using a filtered convective 
velocity in the Navier-Stokes equations |11| . This study has influenced the 
investigation of the Leray-a models of turbulence [12] as well us other mod- 
els of turbulent flows. The Lagrangian Averaged Navier-Stokes- a (LANS-q) 
uses similar filtering to successfully model some turbulent incompressible 
flows [131 El IS] • Although we do not use Lagrangian type averaging in 
this manuscript, it has ultimately inspired the concept of spatial filters on 
equations of compressible flows as a potential model for shock formation and 
turbulence. 

The spatial averaging in the CFB equation has also been extended to both 
the compressible homentropic Euler and full Euler equations in |16| \T7\ IT5] . 
For instance, Bhat and Fetecau [T6] applied a Leray- type averaging with the 
Helmholtz filter g to obtain the equations 

Pt + upx + pux = 

where 

la \ a 

For the isothermal case, 7 = 1, the authors show the averaging prevents the 
crossing of characteristics and prove the unique existence of global smooth 
solutions. They further show that the solutions will converge strongly to 
a weak solution of the homentropic Euler equations as a — )• 0, however, 
they conclude from numerical simulations that this is not the physically rel- 
evant solution (Note: for 7 7^ 1, this Leray-type averaging does not prevent 
characteristics from crossing). 



A promising averaging technique introduced in [17| . called the observable 
divergence method, was applied to both the homentropic Euler and full 
Euler equations in one dimension [TH] . Mohseni ^7\ derived these equations 
from basic principles. Unlike the aforementioned Leray-type averaging, the 
observable divergence method appeared to numerically capture the entropy 
solutions for the shock-tube problem and the Shu-Oscher problem. This 
particular type of averaging will receive full attention within this manuscript 
since it appears to be a valid shock-regularization for the equations of gas 
dynamics. The importance of this regularization is its potential to model 
small scale behavior capable of capturing shocks and (possibly) turbulence 
in a single technique, thus reducing the computational cost demanded by 
many problems in engineering flows. 

Influenced by such work, we shall adopt spatial averaging for the transport 
equations by applying the method of observable divergence to the continuity 
equation and Burgers equation then obtain analogous results to the method 
of vanishing viscosity. This manuscript will be organized as follows. Section 
[2] will review basic material on the transport equations, delta-shock waves, 
and the Riemann problem. Sections [3] and [4] will discuss the class of aver- 
aging kernels in consideration and will introduce and briefly motivate the 
method of observable divergence. A derivation of the observable transport 
equation is provided. Section [5] examines the observable transport equa- 
tions with an existence and uniqueness result in the class of smooth initial 
conditions. The one-dimensional case is considered with a discontinuous 
initial condition consisting of a single decreasing jump since this provides a 
Dirac delta distribution solution for the original transport equations. It is 
shown that the regularized equations share the same Dirac delta distribution 
solution. 



2 Preliminaries 

Let us first review several fundamental ideas and results for the transport 
equations. This section reviews the one-dimensional transport equations 
and its distribution solution for Riemann initial data. 



Pt + {pu)x - 


= 




= 



2.1 The one-dimensional initial value problem 

For now, let us restrict our attention to the one-dimensional initial value 
problem 

(2.1) 

(2.2) 
(p(x, 0), n(x, 0)) = (po(x), no(a;)) (2.3) 

This gives us a simple example of a system of a nonlinear nonstrictly hy- 
perbolic PDE in which classical weak solutions may not exist. Instead, one 
must formulate a broader notion of a distribution solution and include a 
new type of nonlinear hyperbolic wave, a delta-shock wave. To illustrate 
this idea, assume that the initial conditions are smooth and u'q{s) < for 
some point s. Then the characteristic that passes through the point s at the 
initial time t = is given by x = s -|- uo{s)t, and the velocity is constant, 
u = uq{s), along the characteristic. In addition, it easily follows that along 
the characteristic 



and 



This implies that 



9x ' 1-1- UQ{s)t 

Po{s) 



p{x,t) 



l + u'f^{s)t' 



lim {p,u) = {+oo,+oo) 

-«;)(s)-i 



along the characteristic. Hence, we can only define a smooth solution up 
until we approach the finite-blowup time t = —u'f^{s)~^ where the gradient 
and density exhibit singularities. 

2.2 Distribution solutions 

Define a two-dimensional weighted delta function w{s)6c supported on a 
smooth curve C parametrized by t = t{s), x = x{s) {a < s < b) hy 



<w{s)5c,^Pix,t)>= / w{s)tp{x{s),t{s))^/x''^ + t''^ ds 

J a 

for ah V^ G C^(M2_). 

Then we define two distributions in the dual space ((7(j"(M^))* by 

poo f'OO 

<h,(j)>= / / h(pdxdt+ <w{s)Sc,(l)>, (2.4) 

Jo J-oo 



and 

roc foo 



<hu,(j)>= / / hu(()dxdt+ < w{s)6c,(t> >, (2-5) 

Jo J-oo 

for (p G C^(Mi) and h locally integrable. 



Definition 1. Define two distributions G (C^(M^))* as in (2.4) and (2.5). 



Then we say that p is a (Dirac delta) distribution solution to (2.1) if 

< p,(f)t> + < pu, (f)x >= 

for all (p £ Co°(^+)- -^^'" ^^^ Burgers equation, its standard definition of a 
generalized or weak solution will be used. 

2.3 The Riemann problem 



Consider (2.1)-(2.2) with the following initial condition 



(po(x),no(x)) = | ^^^^^^^^ ^^Q (2.6) 

where pi, pr > 0. This is nonstrictly hyperbolic with a repeated eigen- 
value X = u and corresponding eigenvector r = (1,0)^. The eigenvalues 
are linearly degenerate, thus the solutions can be constructed by contact 
discontinuities, a vacuum or delta shock wave connecting the two constant 
states {pi,ui), {pr,Ur). 

case(i) ui < u^ 

A vacuum occurs between two contact discontinuities. A self-similarity 
transform ip = x/t yields the solution 

{{pi,ui), x/t<ui 
(0,x/t), Ui<x/t<Ur (2.7) 

{pr,Ur), x/t > Ur. 

case(ii) u* := ui = Ur 

It easily follows that the left and right constant states are connected by a 
contact discontinuity along x = u*t. 

case(iii) ui > Ur 

This is the case in which a delta-shock formation occurs i.e. the density p 
is the sum of an L°° function and a weighted delta function along the shock 

curve. 

{pi,ui), x/t<a 

{p,u){x,t) = <^ {w{t)6{x — at),us), x/t = a (2.8) 

{pr,Ur), x/t > a, 
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Figure 1: Case ui < Ur 



Figure 2: Case ui > Ur 



where 

a = us = -{ui + Ur), ui> a > Ur, and w{t) = „ i[pu] - cr[p]). 

Note that u is given the value u^ along the shock. This is justified because 
if the values of u are changed on a set of 2-d Lebesgue measure zero, then u 
would still remain a weak solution for Burgers equation. 



3 Filters and averaging kernels 

In this section, we define the acceptable filters used in averaging the chosen 
quantities in the equations. For a given real-valued function /, we introduce 
an averaging kernel g and define the filtering of / by the convolution operator 



f = 9*f 



9{x-y)f{y)dy 



for X belonging to some given domain fi C M". 

We shall only consider averaging kernels g G M^^'^(M") that satisfy the 

following: 



1. 



gd'x. = 1 



2. 



^(x) > for all X € 



|xi| < |X2| =^ \g{^l)\ < |5(X2) 



4. 



|xi| = |X2| => \g{xi)\ = |C/(X2)| 



5. 



lim k^(k) = 0, where g is the Fourier transform of g 

|k|^oo 



In the physical sense, the averaging should only provide non-negative weight 
to particles, have no preferential direction, and should give more weight to 
particles that are physically closer. In addition, we prescribe a parameter, 
a > 0, to such a filter such that 

!,' = -!,(-)■ 

a \a/ 

This parameter a acts as a scaling of the kernel and controls the level of 
filtering. One example of a commonly studied filter is the Helmholtz filter 
f = f — a^A/ corresponding to the averaging kernel 

9{x) = —exp 

2a \ a 



4 Flux Averaging and the Method of Observable 
Divergence 

This section will examine how the particular spatial averaging method is 
applied to the transport equation including background on its development. 

For the Burgers equation, the filtering of the convective velocity i.e. the 
convectively filtered Burgers (CFB) equation ut + uux = 0, provided a valid 
regularization by preventing the crossing of characteristics |10l |8] . With the 
success of this regularization, an extension to this averaging was sought for 
more general conservation laws. One such technique-dubbed the observable 
divergence for a vector field-entails an averaging of the flux term J17j . More 
precisely, if we are given a scalar function / and a vector field V, we define 
the observable divergence of the vector field F = /V by 

odiv{F)=JV-V + V-Vf. (4.1) 

Then conservation laws of the type ft + ^ • (/V) = are modified into its 
'observable' conservation law ft + odiv{F) = 0. Furthermore, Mohseni later 
derived the 'observable' conservation law equations from basic principles 



To explain how the CFB equation inspired this particular averaging, note 

2 

that the fiux term for the one-dimensional Burgers equation, {\)x = {^)x-, 

may be expanded as 

1 

-{uux + uxu). (4.2) 



With this in mind, Norgard and Mohseni [18) developed the observable 
divergence method as a generalization of the averaging used in the CFB 



equation by observing that averaging the non-differentiated terms in (4.2) 

-{uu^ + Uxu) (4.3) 

yields the relationship odiv{^) = 2{uux + Uxu) = uux- Therefore, we may 
formally obtain the CFB equation by applying the method of observable 
divergence to Burgers equation i.e. 

ut + odiv I "^r I = ^^=^ "Wi + ^^x = 0. 

We extend this averaging for Burgers equation in multiple-dimensions by 

Ui -F u • Vu = 0. 



Clearly, (4.1 ) reasonably follows from this concept of expanding the flux and 



filtering the non-differentiated terms as illustrated in (4.2) and (4.3). Now 
consider the continuity equation pt + V ■ (/ou) = and apply the observable 
divergence method to F = (pu) to obtain the equation 

pt + p V • u + u • V/9 = 0. 

We shall call this modified continuity equation the observable continuity 
equation. Hence, the observable transport equations-^the system combining 
the observable continuity equation and CFB equation-is derived by applying 
the method of observable divergence to the transport equations. 



5 Observable equations 

Consider the observable transport equations 

/9t + pV • u + u • Vp = 



ut + u-Vu = 0. ^^'^' 

This section will address the very importannt and desirable property con- 
cerning the existence of solutions to the observable transport equations in 
some given domain O x I C M" x M and satisfying some given initial data 

(p(x,0),u(x,0)) = (po(x),uo(x)). 

The results of this section is summarized as follows. An existence and 
uniqueness result for the multi-dimensional transport equations with smooth 



initial conditions will be established. First, the case for the CFB equation 
is covered. Next, a fixed point theorem is used to show the existence and 
uniqueness of a global broad solution for the observable continuity equation. 
Then it is shown that smooth initial conditions will yield classical solutions. 
The use of this spatial averaging as a legitimate alternative to vanishing 
viscosity requires it to correctly model the delta-shock phenomena. There- 
fore the Riemann problem for the one-dimensional problem is considered 
since it provides the simplest example of discontinuous initial data that pro- 
vide solutions that contains a delta-shock distribution solution. It will be 
shown that the observable transport equations share the exact delta-shock 
distribution solution as the non-regularized transport equations. 

5.1 Existence and Uniqueness 

If we start with smooth initial conditions uo{x) G C^(M"'), we can show 
the existence and uniqueness of the velocity field u{x,t) G C"'^(M",M) to the 
initial value problem 



Ui + u • Vu = 

u(x,0) = uo(x). 



(5.2) 



We briefly describe the procedure to showing this, but the reader is referred 
to |H] for the complete proof. Associate the map (j)t{^) • J^" '~^ J^"' satisfying 

d _ 

— 0f(s) =u((/)t(s),t) 

with initial condition i;^t=o(s) = s. This map can be interpreted as mapping 
the initial p osit ion to its position at time t and is called the Lagrangian 

map. Thus (5.2) yields that — u(</)t(s), t) = 0, which implies that 



u{(ptis),t) = u{(j)t=o{s),0) =uo{s). 

One can show that the Jacobian of (pt, J{4>t), is nonzero for all s £ M" and 
for every t > 0. Hence, the mapping 1 1— )• (/>t(s) is a diffeomorphism for each 
fixed t. From this result, it follows that the velocity field is given by 

u(x,t) = uo(07^(x)). 

Remark: Similarly, if the initial condition uo(x) G L°°{W^), then there exists 
a unique solution u(x,t) G L°°(M",M). 

Now we are only left with determining a solution for the density function, 
p{x, t) to the initial value problem 

pt + pV-u + \i-Vp = , , 

p{x,0) = po{x). ^^-""^ 
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(^, ^l 



x(t;^, xj 



x(0;%, T ) 



Figure 3: Characteristic passing through (5,t) 



Initially, we prove the global existence of a bounded measurable solution that 
satisfies the observable continuity equation along the characteristic curves. 
From there, we show that if we choose the initial condition to be continu- 
ously differentiable and impose further regularity conditions on the velocity, 
then the previous bounded measurable solution turns out to be continuously 
differentiable as well. 



Definition 2. A closed region D = i^ x I C 



xM. is called a domain of 
determinacy for { 5.3) provided that the following holds. For every (^, r) G D, 



the characteristic curve {(x(t;^,r), t)} is entirely contained inside D. 



From above, recall the characteristics are defined to be the map (j)t{s), and 
we showed that these characteristic curves are nonintersecting for different 
initial point s. Moreover, it was shown that the velocity was constant along 
these curves with respect to time. From this point on, we shift our perspec- 
tive slightly. From the implicit function theorem, we can find a C^ mapping 
5 : M ^ M" such that (/>i(s) = Mdit)) =■ ^it) for t > 0. That is the charac- 
teristics may be parametrized on our domain of choice by the curve (x(t), t). 
As a consequence, for every point (^, r) in our domain, we denote x{t] ^, r) 
to be the unique characteristic curve passing through (^, r) (see Figure 3 
below). Before stating and proving our desired result, let us first introduce 
several concepts and definitions. 

So suppose that D is a domain of determinacy, and let X be the Banach 
space of bounded measurable functions v : D >-^ M with respect to the 
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weighted essential supremuin-norm 



-2Lt\ 



\\v\\ = esssupe |t^(a;, t)| 

D 

where L is a fixed bound for V • n on D. 

dx 
Furthermore, the characteristic curves x{t; ^, r) clearly satisfies — = u(x{t),t) 

Thus 

^(x(t; e, T),t) = -p(x(t; e, T),t) V • u{x{t; ^ T),t). 



This observation motivates our definition of a broad solution to (5.3). 



Definition 3. Let H : X ^ X be the map H[p]{x,t) = — p(x, t)V • u{x,t). 
We say an integrable function p : D ^M. is a broad solution for the problem 
(5.3) provided that 

p(C,r) = po(x(0;C,r))+ / ^(x(t;C,r), t,p(x(t; e,r),i)) di 

holds for (^, t) £ D almost everywhere. 

Now the following lemma will prove useful in establishing a Lipschitz-type 
inequality required in Theorem 1 where a fixed point theorem is used to 
show the existence and uniqueness of a global broad solution. 

Lemma 1. Consider again the map H[p]{x,t) as defined previously. Then 

\H[pi]{x,t) - H[p2]{x,t)\ < Le2^*||pi -P2II 
for {x,t) G D and pi,p2 S X. 

Proof. 

e-^'''\H[pi]{x,t) - H[p2]{x,t)\ < Le-^'-'lMx^t) -P2{x,t)\ 

<L esssupe" ^\pi{x,t) —p2{x,t)\ 

D 

< L||pi -P2II < -^ll/Ol -P2||- □ 

Theorem 1. Let D = 0, x L be our domain of determinacy. Suppose that 
po : U ^ M. is a bounded measurable function and uq G W^'°°{Q,). Then 
there exists a unique bounded measurable solution u : D ^ M. to the initial 



value problem (5.3) 



value problem {5.2) and a unique broad solution p : D — )• R to the initial 
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Proof. Define the integral operator T : X — t- X to be 



Tp{^, r) = po{x{0; C, r)) + / H{x{t; C, r), t, p{x{t; C, T),t)) dt 

Jo 

Claim: The map T : X ^ X is a. contraction map. 
For (x, t) G D, pi,p2 G X, and by Lemma 1, 

\Tp,{C,T)-Tp2itr)\< r \H[pi{x{t;tT),t)]- H[p2{x{t-C,r),t)]\dt 

Jo 

f'T f'T 

< / Le'^^^Wpi- p2\\dt< L\\pi- p2\\ e^^^dt 



JO 

< JllPi-P2||(e2^^ 



<IIPi-P2||(^e2^* 



- 2' 



<^e^^nipi-P2| 



It follows that 



That is, 



e ^^^IT/Oi - Tp2\ < -\\pi - P2\\- 

\\Tpi-Tp2\\ < -\\pi - P2\\- 

Hence, by the Banach fixed point theorem [T9], there exists a unique element 
p E X such that Tp = p. The existence and uniqueness of u has been 
discussed previously. D 

Classical Solutions 

Let us now prove that if we have C^ initial data, then the unique bounded 



measurable solutions to (5.1) are in fact within the class of C^ functions. 



Theorem 2. Let po : 0, ^M. and uq : 0, ^M. be continuously dijjerentiable 
and assume D = Q x I is a compact domain of determinacy. Moreover, 
suppose that the second partial derivatives of uq are bounded on (7. Then 



there exist unique solutions u{x,t), p{x,t) € C^{D) to the IVP (5.1). 



Proof. From the previous theorem, p is the uniform limit of the sequence 
p^"^\ for m > 1 defined inductively by 
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or equivalently as the sequence 



Jo 



Since the initial data po is continuously differentiable, we deduce inductively 
that if p^"^' is C^ then p("^+^) is C^ for each n > 0. Since the uniform limit 
of a sequence of continuously differentiable functions with respect to the 
supremum norm is at least continuous, we are only left to show that p has 
continuous partial derivatives on D. 

Take the partial derivative of this with respect to ^i for 1 < i < n to get 



_d_ 



(m+l) 



(e>-) = E 



k=l 






.(V.n)+^V. 



fc=i 



dxk dii 



dt 



where the integrals are evaluated along the curve t i— )■ x(t;^, r). Now for 
each integer m>l and every r > 0, define 



Ym{T) = sup 



dp{^) 



di, 



(e,T) 



and 



■Z'm(r) = sup 






5e* 



d^i 



Then there are constants Ci, C2 such that 



Yrr^+iir) < Ci + C2 / Y„,it)dt. 
Jo 



This inequality implies that 1^(t) < Cie^^. Therefore {1^} is uniformly 
bounded. 

Since T was shown to be a contraction map, we can find a constant C3 such 
that 

|^m+i(^^^) - P™(C,r)| < C32-'" for all {(,t) G D, m > 0. 
With this we have, for some constant C4, the following inequality 

Zm+l < 2^™C3 + C4 Jg Zm{t) dt. 
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This inequality implies that Z^ < T-^iC^ + \)e^^^'' . Hence, Em=o ^™ 

converges uniformly. It follows the partial derivatives -fe^p > q^ uniformly 

therefore q^ is continuous on D. A similar argument can be applied to the 

case concerning -S. As briefly discussed above, the result for u has been 
covered above. D 

5.2 The Riemann Problem 

The existence of solutions to the observable transport equations has been 
covered and it was revealed that the smoothness of the solutions depend on 
the smoothness of the initial data, but it has yet to be verified if the regular- 
ized equations recover the delta-shock solutions, as is normally accomplished 
with vanishing viscosity. The basic Riemann problem is a well-known exam- 
ple of a discontinuous, piecewise constant initial data that provides a Dirac 
delta-shock solution. Remarkably, it will be shown that for arbitrary a > 0, 
the observable equations may contain delta-shock distribution solutions as 
well. In fact, it shares the same distribution solution as the non-regularized 
equations. So consider the Riemann problem for our regularized system 

Pt + upx + pux = 

ut + uux = (5-4) 

{p{x,0),u{x,0)) = {po,uo) 



ipoix),uoi.)) = { [Py:!\ !<° (5.5) 



where 

{pr,Ur), X > 0. 

We shall restrict ourselves to the case ui > u^ since this is the case where 



the solution is a delta-shock wave for (2.1)-(2.2). Recall the characteristic 



curves via the Lagrangian particle map 4't{s). Again, we have 

^=5(0,W.*). (5.6) 

The solutions for the Riemann problem for CFB equation is studied in |20j . 
Here, the CFB Riemann problem is solved using the method of characteris- 
tics instead of using the usual idea of weak solutions involving test functions 
and integration by parts. We cover the necessary results from j20] . It is 
shown that 

Ui+Ur 

where a = . 
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Remark: This solution is independent of a. 



Furthermore, without any restriction on u/ and u^, it foUows that (5.6) 
impUes 

ui + {ur-ui)ip-{^), x/t<a 

u{x,t) = < Us, x/t = a (5-8) 

^ Ur + iUr-Ul)ip+{^), x/t>a. 

Here ^~{x) = J^^ g~{y) dy and ^'^{x) = J^ g^{y) dy where the averaging 
kernel is expressed as a piecewise function 

"W - { :;!:!: : 1 1 (^■») 

Due to conditions placed on the averaging kernel, c^ decay to zero as |x| — > 



oo. Therefore it follows from (5.8) that 



ui, x/t < a 
lim u{x,t) = i us, x/t = a (5.10) 

^^ Ur, X/t > a. 

In the case ui > Ur, this limit recovers the same solution for u{x,t) as that 



in (2.8). 



By the method of characteristics, the correct solution for u is obtained, 
however, we are left to show that the delta-shock solution is preserved. So 
unlike the case for the CFB equation, we must express pt + upx + pux = 
in conservation form in order to define a distribution solution as in Section 
[2j For brevity, let us restrict our attention to the Helmholtz filter, u = 
u — a^Uxx- Choose a test function (f) with compact support and observe that 

Pt + pUx + upx = Pt + {pu + pu)x - Px^ - PUx 

= pt + {pu + pu)x - Px{u - a^Uxx) - (p- a^Pxx)ux 
= pt + {pu + pu)x - {pu)x + a'^{PxUx)x- 

Integrating this against (p and using integration by parts gives us 



p<l)x + pu(l)x + p{u - u)<t)x + a^{pxUx)(px dxdt. 

In light of this observation and in accordance with Section [2j we say that p 
is a distribution solution to the observable continuity equation if 

< p,(i)t> + < pu,(j)x> + < p{u -u),(j)x > + < a^PxUx, 4>x >= 0. 
The next theorem will show that the original delta-shock distribution solu- 



tion to the transport equation also satisfies the regularized equations (5.4) 
as well. 
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Theorem 3. Let a > and suppose that ui > u^, a = "'^"'' E {ur,ui), and 
[pu] —cr[p] > 0. Furthermore, let g be the Helmholtz filter. Then {p,u) given 
by 

p = h + w{s)5c, (5.11) 



where 



h{x,t) 



and 



C := {t = s,x = as : < s < oo}, 



ti/, xjt < a 
u{x,t) = ^ Us, x/t = a 

Ur, x/t > (T, 



(5.12) 



(5.13) 



satisfy the Riemann problem to (5.4)- 



Proof. Recall that u(x, t) and u{x, t) was determined using the method of 
characteristics. For the Helmholtz filter, it follows that 



p(x, t) = h + w{s)5c < 



Pi - \{pi - Pr)exp{^^-^), x/t < a 

w{t)6{x — at), x/t = a 

Pr -\{Pr- pi)exp{-^^^^), x/t > a, 



(5.14) 



u{x,t) 



ui - ^{ui - Ur)exp{^^-^), x/t < a 

Us, x/t = a 

Ur -\{ur-ui)exp{-^^^^P^), x/t>a, 



Px{x,t) 



2a 



h(Pi - Pr)exp{^^-^), x/t < a 
), x/t > a. 



(Pr - Pi)exp{- 



(x—ai) 



and 



-i^{ui - Ur)exp{'^-^), x/t < a 
^{ur - ui)exp{-^^^^), x/t > a. 



(5.15) 



(5.16) 



(5.17) 



Ux{x,t) = 

To complete the proof, it is left to show the following claim holds true. 
Claim: < p,(f>t > + < pu, (px > + < p{u — "u), ^x > + < o.'^Px'^x, 4>x >= 
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To prove this claim we compute each term on the right-hand side then verify 
that their sum is zero. 

(i) < p, 0t > + < pu, <px >= / / p(t>t + {pu)(t)x dxdt 

oo /'oo /'oo px/a 

/ pi(t)tdtdx+ j j prCptdtdx 

Jx/a Jo Jo 

POD 

+ / w{t) ((/)f (cri, t) + us(px{(^t, t)) \J\ + cr2 dt 

Jo 
+ / / pi{ui- -{ui - Ur)exp ( j 4>x dxdt 

r r , 1, , [-{x-at)\ ^ ^ ^^ 

+ I / Pr(Ur {Ur — unexp (px uxdi 

Jo Jat 2 \ a J 

""^ '■"* 1 / ^ fx-at\ ^ ^ ^ 

-pi{ui - Ur)exp cpx dxdt 

10 J-oo 2 \ a J 

""^ '""*!, , f-{x-at)\ ^ , , 

-pr[Ur - uijexp (px dxdt. 

10 J-oo 2 \ a J 

Similar calculations for the other terms gives us 

_ _ f^ Z"'^* 1 / X — at\ 
(u) <p{u-u),(j)x>= / / -piiui - Ur)exp I ](f>xdxdt 

r°^ ro-t 1 /2(x - at)\ 

-^{Pi - Pr){ui - Ur)exp (px dxdt 

10 J-oo 4 \ a J 

f^ f^ 1 ^ , f x-at\ ^ ^ ^ 

+ / / 7,Pr{Ur - Ul)exp (px dxdt 

Jo Jat 2 V " / 

°^ f°^l, ,, , { 2{x-at)\ , , , 
/ -{pr - Pi){Ur - ui)exp \ \(pxdxdt 

Jat ^ V « / 



and 

2- - , f^ r* 1/ ^/ N /^2(x-at)\ , , , 
(m) < a p^Ux, (px>= I j jKPi - Pr){ui - Urjexp ( I (px dxdt 

""^ f°^l, ,, , / 2(x-at)\ , , , 
/ -APt - Pi)[ur - ui)exp \ \(pxdxdt. 

Hence, summing up these terms gives us zero i.e. (i) + (ii) + (iii) = D 

For completeness, let us briefly discuss the case for ui < Ur- Recall that the 
original solution is a rarefaction wave with a vacuum and does not contain a 
delta function. The CFB equation will capture a non-entropic weak solution 
of Burgers rather than the rarefaction wave. This is remedied by replacing 
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the Riemann initial condition with a filtered approximation. That is, intro- 
duce filtering to the initial condition independent from the CFB equation, 
then the CFB will recover the entropy solution for that filtered initial con- 
dition [20]. On the other hand, Norgard and Mohseni [^ showed the CFB 
equation recovers the entropy solution for a class of smooth initial conditions 
and proposed a more general program in recovering the entropy solution for 
discontinuous initial conditions. Such a program can prove useful in numer- 
ical simulations. More specifically, by spatially averaging Burgers equation 
via the CFB equation in conjunction with the discontinuous initial condi- 
tion, the entropy solution is captured for the discontinuous initial data as 
filtering vanishes. As a caveat, this follows from |9j if one assumes that the 
CFB equation recovers the entropy solution for arbitrary continuous initial 
conditions. This behavior indicates how the initial smoothness can influ- 
ence how effective spatial averaging will model our equations and suggests 
that the method of observable divergence may be more appropriate for more 
regular initial data. 



6 Concluding Remarks 

A spatial averaging technique-previously applied to the Burgers and Eu- 
ler equations as a valid shock-regularization-was applied to the transport 
equations. Theory developed for the CFB equation combined with the past 
success of numerical analysis for the regularized equations was the primary 
motivator in studying the effects of averaging the transport equations. In 
particular, the method of observable divergence was introduced and applied 
to the transport equations. It was shown that utilizing this carefully chosen 
averaging technique will exhibit a unique classical solution for smooth ini- 
tial conditions. But unlike the method of vanishing viscosity, the observable 
equations will preserve the original delta-shock solution for the Riemann 
problem regardless of the amount of filtering used. Nonetheless, the delta- 
shock is captured when the level of filtering vanishes. These results reinforce 
the validity that this technique is a legitimate alternative to the method of 
vanishing viscosity-further establishing its capability as an effective tech- 
nique in properly modeling general conservation laws. 
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